Entanglement enhancement for two spins assisted by two phase 

kicks 



Gennadiy Burlak^, Isabel Sainz^, and Andrei B. Klimov^* 

^Centra de Investigacidn en Ingeniena y Ciencias Aplicadas, 

Universidad Autonoma del Estado de Morelos, Cuernavaca, Morelos, Mexico 

O ■ Departamento de Fisica, Universidad de Guadalajara, 

O , 

CN ! Revolucion 1500, Guadalajara, Jalisco, 444^0, Mexico. 

C . 

^ . Abstract 

I We study the entanglement dynamics in a two-spin system governed by a bilinear Hamiltonian 

and assisted by phase kicks. It is found that the application of instant kicks to both spins at 
Q_i' some specific moments leads to enhancement of entanglement. This procedure also improves the 

^ , transient character of entanglement leading, for large spins, to a formation of a plateau for the 

a _ 

^ I /-concurrence. We have numerically investigated the spin-spin dynamics for several values of spins 

and observed a substantial enhancement of entanglement in comparison to the evolution without 
^ . kicks. 
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I. INTRODUCTION 



Recently, a significant amount of tlieoretical and experimental research has been carried 
out for generating, manipulating and preserving entanglement, which is the key ingredient 
for developing quantum technologies. The main attention is usually paid to studying en- 
tanglement properties of qubits (two-level systems) (see e.g. 

113) 

, considered basic 

physical resource. Although, higher dimensional systems (qudits) could be more appropriate 
for certain physical applications, entanglement dynamics, generation and manipulation in 
higher dimensions is not as well studied as in two- dimensions. 

In particular, it was suggested in that it is possible to increase the security, the bit 
transmission rate, or both in quantum key distribution protocols by increasing the dimen- 
sionality of the involved systems. The possibility to perform quantum computation based 
on qudits cluster states {^(has also been noted) . Nevertheless, just few years ago some 
studies on entanglement dynamics in higher dimensions have been done. For example, the 
dynamics of a initially entangled two-qudit state was studied in js], while in Q it is shown 
that two non-interacting, initially separable spins can get entangled via a common purely 
dephasing environment. More general results about environment mediated generation of 
entanglement in higher dimensions can be found in [3]. Other aspects of entan^ement in 
higher dimensional systems, like entanglement concentration for two qudits in |8j], and its 
experimental implementation for two qutrits jo] were reported. It was also shown in [lo| that 
the quality of entanglement transfer in spin chains actually increases with the dimension of 
the spin. 

On the other hand, higher dimensional systems (real, as high excited nuclear spins or 
effective, as Dicke-like states) are of interest by themselves, especially in the context of 
studying the quantum-classical transition, when the system's dimension becomes sufficiently 
large. However, it is well-know that there are no reasons to expect that the classical behavior 
can be approached in the large dimension limit for every state HI]. In particular, the 
coherent spin states (CSS), IC'S'), jl2| are the most classical ones, and some specific linear 
combination of them ma y re veal non-classical features, like entanglement, even in the limit 
of large values of S (see [l3|). 

One of the regular ways for generating and measuring entanglement in many spin 1/2 
system is through spin squeezing (see pLJ] and references therein), where quantum correla- 
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tions naturally appear in the basic concept of squeezed spin states (SSS) [15]. The simplest 
Hamiltonian which produces squeezing starting with coherent states located on the equator 



of the Bloch sphere is H = gSl, where = X]j '^i^'V^ is a collective spin operator 15|. 

Unfortunately, there are two problems with such dynamical generation of squeezing and 
entanglement: a) the squeezing (and entanglement) is always transient, i.e. it appears 
periodically ; b) it does not reach its maximal value, except the simplest two-qubit case. 

In the case of large two spins the situation is rather similar. It is possible to generate in an 
easy way the entanglement between these spins, but it would resemble the above mentioned 
disadvantages as in the symmetric combination of many qubits. 

In the present article we address the following question: is it possible to improve the 
entanglement properties of a two qudit system in a relatively simple way? 

We will show that by applying instant kicks to both spins at some particular moments 
we can not only enhance the entanglement, but, more importantly, essentially improve its 
transient nature, i.e. reduce the distance between the maximum and the minimum values, so 
that the system evolves into an approximate "steady" state of the interaction Hamiltonian 



arge spins. We also discuss its possible 



with a high value of entanglement for sufficiently 
application to the purely dephasing environment jo. 



II. THE MODEL 

Let us consider two spins 5*1 and 5*2, of dimensions 25*1 + 1 and 25*2 + 1 respectively, which 
interact according to the Ising-like Hamiltonian, 

H = gS,iS,2. (1) 

This kind of interaction was proposed to implement one way quantum computation with 
many-level cluster states (see J] and references therein). Suppose that each spin is initially 
in the eigenstate of Sxi with zero phase, S^i \Ci = 1, Si) = Si IQ = 1, Si), which is a CSS 
placed on the equator of the Bloch sphere, 

Ki — — Oi " 

where i = 1,2 and Szi \ki. Si) = ki \ki. Si), so that the joined state for the system composed 
by spins 5*1 and 5*2 is |^o) = iCi = 1, 5*1) IC2 = 1, 5*2) . 



Since |\l/o) is not an eigenstate of the Hamiltonian, in the course of evolution the joined 
state 

= U{t) l^o) , U{t) = e-**»^^l^^^ (2) 
becomes entangled at some instants. As an entanglement measure we will use the (normal- 

n 

ized) /-concurrence [16|, that for pure states takes the form, 

'^^ = ^;^(1-Trpi(2)), 

where pi(2) = Tr2(i)(p(t)) is the density matrix of the composed system p{t) = 

traced by the spin 5*2 (>S'i), and S = min(S'i, 52), such that < C/ < 1, being zero for 

separable states and one for maximally entangled states. 

For the state (E]), C/ = {2S + 1)(1 - V)/2S, where the purity V = TrpUt) is 



V 
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[{2S^)\Ycos''^{gt{k-l)/2) 



From now on we will focus only on the symmetric case: Si = S2 = S. Let us consider 
the two simplest cases: S = 1/2, 1 . For two qubits, S = 1/2, Cj = siv? gt/2, i.e. it oscillates 
between zero and one, so that the spins are maximally entangled at times gt = (2n + l)7r and 
they are separable at times gt = 2mT, n = 0,1, . . .. Nevertheless, for 5* = 1, the situation is 
already different, since C/ = 3 (4 — 2 cos gt — cos^ gt — cos^ gt) /8, reaches its maximum value 
^max ^ g gg timcs gt ^ 2.2,4.1 + 2mT, and thus the spins will never become maximally 
entangled, although they are disentangled at 5ft = 2?T,7r, n = 0, 1, . . .. The entanglement 
dynamics have a quasi-periodical behavior, as is shown by the (green) dash-dotted line in 
Fig. [11(a). 

In the large spin limit, S" 1, we can approximate the time average of Cj as follows, 

iC, . !^ f ^ - (3) 



2S VV2^ 27r5. 

giving (C/) ~ 0.54 for S = 1, which is quite close to the exact numerical value (C/) = 0.58, 
even if the approximation is done for 5* ^ 1. 



III. ENHANCING ENTANGLEMENT 



In order to enhance the maximum achievable entanglement let us apply instant pulses 
(kicks) to each spin at certain times tk- Such kicks correspond to rotations around the y-axis, 
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(a) 



- S=1 0, with phase kick 
' S^IO, no phase kick 



(b) 



FIG. 1: (Color online) Time evolution of Cj, without kicks (green, dash-dotted line) and with 
phase-kicks (blue, continuous line) for (a) S = 1, with the numerical optimization times ti = 1.6 
and t2 = 3.9, and (b) S = 10, in this case ti = 1.9 and t2 = 4.2 where the inset shows the details 
for the second phase kick. 

and for the j-th spin are represented by the operators Rj = e"*'^'^"/^. The main idea is as 
follows: initially both spins are in eigenstates of S^i, which can be represented by localized 
distributions on the equator of the corresponding Bloch spheres. In the course of evolution 
governed by Eq.([l]) the state become partially entangled. At some appropriate moment ti 
we apply a 7r/2-rotation to both spins so that each component of the angular momentum 
basis \ki, S) instantly transforms into an eigenstate of S^i- Afterwards, the evolution under 
the Hamiltonian ([1]) is continued for a while, which continues coupling both spins. Then, at 
some time ^2 we apply the inverse it/ 2 -rotation that transforms again each element \ki, S) 
into an eigenstate of Sxi, after that we allow the system to evolve further: 

\^{t)) = U{t - t2)R-'U{t2 - h)RU{h - to) l^o) , 

where R = Ri ® R2. The instants ti and t2 are chosen in order to optimize the maximum 
value of achievable entanglement after the second kick. A similar procedure was proposed 

n 

to optimall y cr eate squeezing in Dicke states [17J and entanglement for continuous-variable 
systems in 
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19|. 



The reason for applying 7r/2-rotations can be explained as follows: we want to correlate 
as many components of the spin states as possible, which cannot be done only with the 
Hamiltonian evolution, since the number of components of the angular momentum basis, 
initially involved in the CSS, is of order v^- By applying the kicks at some specific moments 
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FIG. 2: (Color online). Distribution of the two-spin probability as function of the 

number of state k in the long basis {0, . . . , (25* + 1)(2S' + 1)} for S = 10. (a) for the initial state, 
t = 0, (b) before the first phase kick at t < ti, (c) after the first kick, but before the second kick 
at ti < t < t2, and (d) after the second kick t > t2- 



we are able to spread the distribution making it more uniform |18l.ll9[|. while the Hamiltonian 
evolution produces correlation between spin components in the angular momentum basis, 
leading to entanglement. 

In Fig. [21 we plot the distribution of the two-spin probability (\E'(t)|\&(t)) as a function 
of the state number k, corresponding to the basis \ki,10) |A;2,10), in the following way, 
| — 10, 10) | — 10, 10) corresponds to = 0, the next state |— 10, 10) |— 9, 10) corresponds to k = 
1, and so on. For the initial state (Fig. [2] (a)), the distribution is centered on |0, 10) |0, 10), 
that is, k = 2S{S + 1) + 1 = 221. The Hamiltonian evolution does not affect the distribution 
of the amplitudes of the components of the basis k {\ki, S) \k2, S)) as shown by Fig. [2] (b). 
It can be seen in Fig. [2] (c) that after the first instantaneous rotation of each spin, the 
distribution of (\E'(t)|\E'(t)) spreads over the states k, this effect is increased after the second 
instantaneous rotation of each spin (Fig. [2] (d)). 

For semi-integer spins (5* = {2n + l)/2, n = 0,1, . . .), the situation is slightly different: 
the initial distribution is similar to Fig. [2] (a), but instead of having a single maximum at 
the center, it has two maxima, which leads to non-essential (but observable) differences in 
the dynamics. 

The analytical expressions are very cumbersome for higher spins, so that we have to use 
numerical optimization for 5 > 3/2. 
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FIG. 3: (Color online) The behavior of (C/) against the spin number S < 12. The (green) squares 
joined by the dashed line correspond to the two-spin system under the Hamiltonian evolution given 
in Eq. [TJ The (blue) circles joined by a continuous line correspond to the enhanced behavior, for 
times t > t2. 

As noted before, for the simplest case, when S = 1/2, the system oscillates between the 
separable and maximally entangled states and the kicks produce just a phase shift, that 
does not affect the dynamics. 

For S = 1, we have been able to analytically optimize Cj, and found that the entanglement 
acquires its maximum value when the kicks are applied at times ti ^ 1.6 and ^2 ~ 3.9 (which 
is also in accordance with the numerical results). For these times, Cj reaches the maximum 
^max ^ g gg times t ~ 7.1 + 27m, which should be compared with that obtained by the 
pure Hamiltonian evolution, where Cf^^ ^ 0.88. However, it is more important that after 
the second kick the quasi-periodical behavior of Cj has the minimal value C™™ ^ 0.57, which 
means that the spins are entangled at any time t > t2 > ti, a.s shown in Fig. [1] (a), and (Cj) 
increases from (C/) ~ 0.58 to (C/) ~ 0.71. 

In Fig. [T] (b), we plot the entanglement dynamics for S = 10. In this case, the effect of 
enhancing the maximum of entanglement for t > t2, is not that pronounced as for = 1, it 
increases from Cf^^ ^ 0.90 to C™^^ ^ 0.96. However, another important effect emerges: the 
transient nature of entanglement becomes significantly improved, i.e. after the second kick 
Cj rapidly oscillates with a very small amplitude around its average value (Cj) ~ 0.95. So 
that the spins are entangled for all times t > t2, with a high value. On the other hand, we 
observe, that in the case of pure Hamiltonian evolution Cj periodically oscillates between its 
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maximum Cf 



0.90 and zero. 



The average Cj grows slowly as S increases, and even if the maximum and minimum 
are getting closer, the rate of increase of the average is not as fast as could be expected. 
This is a consequence of the growing character of (Cj) for the Hamiltonian dynamics for 
large S as seen from Eq. To represent such behavior clearly we plot in Fig. (jS]) (Cj) 
for the Hamiltonian evolution (squares) and after applying the sequence of kicks (circles), 
numerically obtained. 

IV. DISCUSSION 

We have studied the dynamics of a two-spin entanglement assisted by two phase kicks. 
The instantaneous kicks were applied at times such that the maximum of Cj has been 
(numerically) optimized. We found that such kicks not only lead to an enhancement of the 
two-spins maximum entanglement, but also to a substantial increasing of the minimum and 
average values of Cj for times after the second kick. Such an effect becomes much more 
pronounced for large spins S" 1, when C/ tends to a constant value very close to unity. 

It is worth noting that the interaction Hamiltonian similar to ([1]) appears in sev- 
eral effective processes. In particular, when two non-interacting spins are embedded in 
a dephasing environment \(y\. The Hamiltonian describing such evolution has the form 
H = QiSzi + ^2Sz2 + Y.k^ka\ak + {S^i + Sz2)J2k9k{ak + 4), such that Uk > ^1,2, 
Qk- The fast field is effectively decoupled from the spin subsystem and thus can be 
adiabatically eliminated in a standard way [2^ by applying a small unitary transforma- 
tion U = exp[{Szi + Sz2)J2k^kiO'k ~ (^k)], Cfc = Qk/^k, SO that the effective Hamiltonian, 
ife// = UHU'^ (diagonal on the field variables) takes the form. 



where e = 'Ylikdk^k- It can be observed that the effect of bosonic bath in this case is 
reduced to an effective interaction between spins, which may give rise to a creation of spin 
entanglement and a quadratic phase shift that is not important for entanglement generation. 




{VLi - eSzi)Szi + (SI2 - eSz2)S: 




k 
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